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Abstract 


The design, alignment, calibration, and field deployment of a solid-state lightning detector is 
described. The primary sensing component of the detector is a potassium dihydrogen phosphate 
(KDP) electro-optic crystal that is attached in series to a flat plate aluminum antenna; the 
antenna is exposed to the ambient thundercloud electric field. A semiconductor laser diode (A = 
685 nm), polarizing optics, and the crystal are arranged in a Pockels cell configuration. 
Lightning-caused electric field changes are related to small changes in the transmission of laser 
light through the optical cell. Several hundred lightning electric field change excursions were 
recorded during five thunderstorms that occurred in the summer of 1998 at the NASA Marshall 
Space Flight Center (MSFC) in northern Alabama. 
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1. Introduction 


The most common means for measuring lightning-caused changes in thunderstorm electric 
fields involves the application of a flat metal plate antenna and a charge integrating apparatus. 
Because the normal field E n at the surface of an exposed antenna plate satisfies the normal 
component boundary condition E n = QI(As o ), where Q is the induced charge, A is the plate 
area, and £ 0 is the permittivity of free space, a change in the thunderstorm field must be 
accompanied by a change in induced charge. The integrating apparatus used to measure the 
changes in induced charge normally consists of an electronic integrating circuit composed of an 
electrometer and having in the feedback loop a resistor R and capacitor C in parallel. The 
inverting input of the electrometer is attached to the plate and the noninverting input is grounded. 
Typical RC decay times of the electrometer output are measured in seconds (slow antennas ) and 
from milliseconds to microseconds (fast antennas). 

From the point of view of the experimenter, the RC decay time should be made larger than 
any time variation of interest [1]. For instance, if one is interested in lightning return stroke 
fields, the RC should be in milliseconds and if interested in total lightning field changes, the RC 
should be in seconds. The decay assures that the sensor output is nullified before another event 
occurs; triggering from the “zero line” is an advantage because one can devote the full vertical 
amplitude resolution and dynamic range of the sensor to the event. The systems described by [2] 
employed RC values of 4 s, 0.3 ms, and 70 /js to look at various temporal regimes in the 
lightning signature. 

In the past 20 years, several new scientific and commercial applications of electro-optic 
technologies have arisen for measuring electric voltages and fields ([3] - [10]), and some of these 
investigations were specifically directed toward the study of thunderstorm electric fields. 
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Electro-optic crystals have offered a convenient and affordable means of controlling the 
intensity or phase of optical radiation in an ever-expanding number of applications that include 
high-speed digital communication, optical beam deflection, and spectral tunable filters [ 11 ]. 
Because the crystals respond quickly to applied voltages and no moving parts are required, 
modulation frequencies of several gigahertz are possible. Liquid crystals have also been used in 
high speed optical shutter applications, and [ 12 ] has used liquid crystals to measure lightning 
electric field changes. 

In the present study, we describe an electro-optic lightning field change detector that was 
recently developed and deployed in the summer of 1998 at the NASA Marshall Space Flight 
Center (MSFC) in Huntsville, Alabama. The sensor is an engineering model composed of both 
fabricated and vendor-supplied components, and the optics were mounted on a small laboratory 
optical breadboard. The RC decay was adjustable from 1 s to 1 fjs and was typically set to a few 
tens of milliseconds during deployment in order to resolve individual strokes in a lightning flash. 
The main sensing components are a flat aluminum plate and a long cylindrical, high impedance 
electro-optic crystal of potassium dihydrogen phosphate (chemically written as KH 2 PO 4 , or 
“KDP” for brevity); some dielectric and optical properties of KDP are provided in [13]. The 
details of the design, alignment, calibration, deployment, and data results are discussed below. 
We begin our discussion with a review of the linear electro-optic effect since it is fundamental to 
the design characteristics of the lightning detector. 

2. The Linear Electro-Optic Effect 

Friedrich Pockels (1865-1913) was the first to describe the linear electro-optic effect in two 
early publications ([14], [15]). Figure 1 shows a typical Pockels cell arrangement for an electro- 
optic modulator. The central component is an electro-optic crystal, KDP. Light from a diode 
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laser (input beam) passes through two polarizers, a quarter wave retarder, and the crystal before 
exiting the optical system. The output beam would typically be measured with a photodetector 
and transimpedance current amplifier system. 

Equations for the transmission of light through an anisotropic medium such as the KDP 
crystal may be derived from Maxwell’s equations and energy conservation. Conductivity cr is 
negligible for the materials we are using at optic frequencies. The electric displacement D and 
magnetic induction B are related through [16] 

V x F = -dfl/dt , VxH = <?D/<?r + J, V-D = p, V-B = 0, (1) 

where F and H are the electric and magnetic field intensities, respectively. [Note: we have 
used “ F ” to denote the laser light electric field here rather than the customary notation “ E ” to 
avoid confusion, i.e., “ E ” and its components denote the (enhanced) atmospheric field at the 
sensor antenna plate as suggested in section 1 above and in the discussion that follows.] The 
vector J is current density, p denotes charge density, and we have the relations D = £ ■ F and 
B = //H. The electric permittivity £ is a tensor for anisotropic materials; it may be brought 
into diagonal form by a suitable rotation of coordinates. We will approximate p by p n , the 
magnetic permeability of a vacuum, for all materials in the optical circuit. 

Energy conservation is described by 

|^n-(FxH) = -^|c?v(F-D + H-B)/2- [</v F J, (2) 

where n is the outward normal to a surface S which bounds a volume v . We will assume that 
p and J are both zero. The equation for energy conservation can be used to show that the 
permittivity tensor is symmetric. 

Suitable manipulations of ( 1 ) for isotropic media yield a wave equation 

(v 2 + & 2 )f = 0. (3) 
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The same equation is satisfied by H. Plane wave solutions pe ,kr are admitted, where p is a 

polarization, k = A:k, k = 2nj A, and k is the direction of propagation. Angular frequency is 
denoted by co , t is time, r is a position vector, and A is the wavelength. With such assumed plane 
wave forms for F and H in anisotropic media, we have the relations H = (k x F)/(rOyt/) and 
D = — (k x H)/<y from which we may obtain (see [17]) 

(k • F)k -k 2 F = -cd 2 ^is ■ F. (4) 

We let k = con/c, where n is the wave index of refraction and c is the speed of light, and refer to 

a coordinate system in which e is diagonal. The form in (4) is that of an eigenvalue equation. In 
general, there will be two nontrivial solutions. 

The following coordinate system will be used. The Z-coordinate will be along the path of 
the collimated laser beam with the direction from the laser to the photodetector considered as 
positive. The zenith direction vertical to the optics breadboard, on which the components are 
mounted, will define T, and X is chosen to complete a right-handed system. Furthermore, we 
define a system (x, y, z) that is identical to the ( X , Y, Z) system but is attached to the crystal; the 
(x,y, z) system defines the crystal principal axes when no voltage is placed across the crystal. In 
the present discussion of perfectly aligned optics, there will be no distinction between the two 
coordinate systems. Later however, in our discussion of misaligned optics provided in section 4, 
the (jc, y, z) system, or “crystal system”, will be oriented differently than the fixed (X, Y, Z) 
“breadboard system”. 

The indices of refraction and their associated allowable polarizations may be obtained 
geometrically. By considering surfaces of constant energy density, an index ellipsoid 

x 2 /« ( 2 +r/V +z 7 n z = 1 ( 5 ) 
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may be constructed [ 11 ]; n x , n y , and n, are the refractive indices along the coordinate axes. 

The allowable polarizations and indices of refraction may be obtained by finding the axes of the 
ellipse formed by the intersection of the index ellipsoid and a plane normal to the direction of 
propagation. 

The index ellipsoid of the crystal may be modified by applying an electric field to the 
crystal. The general formula is given by [1 1] 

(»/ 2 +*i +(v 2 +r 2k e k)y 2 +(«/ 2 +'-3 (6) 

+ 2 r Ak e k yz + 2 r 5k e k xz + 2 r 6k e k xy = 1, 

where the summation convention has been employed. The e k are components of the applied 
electric field, and the r jk are the electro-optic coefficients of the crystal. It is intended that the 
modulating field be along z; only r 63 is needed since r j3 = 0, 1 < 7 < 5 for KDP. 

We assume that the crystal is aligned so that its axes of symmetry are parallel to the 
coordinate axes. The KDP crystal is uniaxial which means that n x =n y = n 0 , where n () has a 
numerical value of 1.507 and is called the ordinary index of refraction. The extraordinary index 
of refraction is defined by n e = n. = 1.467. Since n e < n 0 the KDP is additionally specified as 
being a negative uniaxial crystal. 

When a field e z is applied along the z-direction, the allowed polarizations are 

(-x + y)/V2 and (x + y)/V2. Light polarized along y upon entering the crystal is split into 
these polarizations, the waves along which propagate at different speeds. For a crystal of length d 
having a voltage V = e,d , the phase lag is 

r = - Hj V s X n " r * v ’ (7) 
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where n s (slow) and n f (fast) are the indices of refraction along the allowed polarizations. 

Waves traveling along the z-axis with no modulating voltage have a speed corresponding to the 
ordinary index. Waves traveling in the xy plane with polarization in the z-direction have a speed 
determined by the extraordinary index of the crystal. The crystal’s utility lies in its ability to 
introduce a variable phase retardation between the two waves in response to an applied voltage. 
The linear relationship provided in (7) is the Pockels effect [18]. 

If the laser light field immediately following the first polarizer isFJ, coscot Y, then the field 
of the output beam will be 

F = F 0 [cos(cu/ - r) + cosryf]X/2. (8) 

Here, the phase change due to the path length between optical components is suppressed. 
Reflection losses at the interfaces are relatively small, and are customarily neglected. Exact 
solutions that account for such losses are available ([19], [20]). 

Using (8) to obtain the square magnitude F 2 , the power is obtained by taking the time 
average of F 2 [21]. Dividing by the input power, which is proportional to F Q 2 / 2, yields the 
optical transmission T given by 

T = sm 2 T/2 = sm 2 [(trV)/(2V K )], (9) 

where V„ =2./(2« 2 r 63 ) is called the half-wave voltage, i.e., V = V n yields a retardation of n 
which is one half of 2 n, or “one half wave”. 

Note that the derivative, T'(V) = dT/dV , or electro-optic sensitivity, is zero when there is 
no voltage across the crystal, i.e., r'(0)=0. In order to have the largest possible change in 
transmission with applied voltage for small voltages, it is convenient to start at a point where 
T'( 0) is as large as possible. This may be achieved by introducing an additional phase 
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retardation of /r/2 by inserting a quarter wave plate before the second polarizer. In such a case, 
the transmission of the system is given by 

r = sin 2 [r/2 + ^/4]. (10) 

This process is the optical analog of biasing an electronic component. 

The preceding discussion has summarized the basics of the linear electro-optic effect in 
anisotropic media. Some early experimental investigations of this effect are provided in [22]. 
Note from (9) and ( 1 0) that a high frequency modulation voltage can be applied to the crystal to 
vary the optical transmission and hence the optical amplitude of the output beam. In this sense 
the Pockels cell configuration can serve as a high speed optical shutter. 

However, the Pockels cell can also be used in the reverse sense, that is, a high gain 
measurement of small changes in the optical amplitude of the output beam can be used to infer 
the changes in crystal voltage. This is the operational principle of the electro-optic lightning 
detector developed in this work. The lightning places a transient voltage across the electro-optic 
crystal thereby changing the optical transmission of laser light through the cell. Specific details 
of the sensor follow. 

3. Sensor Description 

Figure 2 provides a simplistic overview of the electro-optic lightning detector developed at 
NASA MSFC. The detector is placed outside on a flat surface sufficiently far from tall grass or 
pointed objects that could produce local field-limiting corona discharges. Note that a 
thunderstorm separates positive and negative charges and thereby generates strong electric fields 
in the atmosphere and on the ground. When a cloud-to-ground flash occurs, as shown in Figure 
2, tens of coulombs of negative charge from the lower regions of the thundercloud (i.e., ~7 km 
altitude) are deposited to Earth. The large charge transfer is accompanied by a substantial change 
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in the electric field A E at the ground, i.e., a A E g > 5 kV/m is not uncommon for discharges 

that are within a range of a few kilometers. Those lightning events that do not transfer charge to 
ground but transport charge wholly within the atmosphere, called cloud flashes, also produce 
substantial values of A R . In general, the field change at the height of the antenna plate differs 
from that at the ground, and the presence of the sensor itself enhances the local field change. 
Figure 2 indicates the net field change at the antenna plate as A E , and this is the stimulus that 
produces a change in crystal voltage given by AV = A e z d . 

When the sensor detects the lightning-caused field change A E, the photocurrent output 
from the optical system is converted to an analog voltage using a Stanford Research Systems, 
Inc. current preamplifier (model SR570), and the analog voltage is passed to a Pentium computer 
where the signal is digitized at 125 Hz using a Data Translation A/D board (model DT2833). 

A closer look at the sensor’s optical system is provided in Figure 3(a); it is based on the 
linear Pockels cell configuration discussed in section 2. The totality of optics are mounted on an 
optical breadboard that is 30.48 cm by 60.96 cm in dimension. The intent was to build a 
reasonably small engineering model that still provided “on-line” adjustability and 
interchangeability of optical components. Future designs of the sensor would be scaled down 
considerably (i.e., the laser and all optical elements could be precision mounted within a single 
sealed tube only a few centimeters long and wide). 

Working from left to right in Figure 3(a), the first element is a 30 mW semiconductor laser 
diode that has a wavelength of 685 nm and an elliptic beam cross section (major and minor axes 
of 4 mm and 1 mm, respectively). We have also used a 4 mW, 670 nm laser diode, but at the 
expense of reduced sensor sensitivity. The beam height above the optical breadboard is 75 mm 
and the beam divergence is negligible across the length of the breadboard. The only 
considerations during operation are to allow the laser a few minutes of “on time” prior to taking 
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data to assure power stability, prevent reflections back into the laser aperture, and prevent 
excessive temperature changes in the air surrounding the laser. 

The next (optional) optic is a beamsplitter that sends a portion of the laser light to a 
photodetector (shown in the background). These two components allow one to monitor the laser 
stability, but because our Melles Griot, Inc. lasers had minimal (< 1%) power fluctuation, the 
monitoring system was not a strict requirement. 

The next element is a post-mounted, vertically oriented Melles Griot dichroic sheet 
polarizer followed by the KDP crystal. The polarizer has a clear aperature of 26.5 mm and a 
thickness of 2.5 mm. The hygroscopic and water soluble KDP crystal is housed in a dry nitrogen- 
flushed air- tight Plexiglass container that is mounted on a two-axis holder; the holder has 
alignment micrometer adjustments in both yaw and pitch. A durable strap fastens the 
crystal/container to the holder and the holder is bolted to the optical breadboard. The crystal has 
a cylindrical geometry (40 mm length, 20 mm diameter) with ring electrodes wrapped around the 
outer crystal diameter at each end. The ring electrodes allow for a larger crystal aperture than 
commonly used washer-type electrodes that are placed directly on the crystal entrance/exit faces. 
Additional details and considerations of ring electrodes are provided in [23]. Direct 
measurements of crystal impedance and capacitance give the values R = 194 GQ and C = 5.51 
pF, respectively, or a crystal RC decay time of ~1 s. (In addition, we have also experimented 
with a post-mounted 3 mm thick, 25.4 mm aperture crystal of potassium didueterium phosphate 
(KD*P) similar to the crystal used in [24] and [25], This crystal had transparent indium tin oxide 
(ITO) electrodes, and, because of its thin width it had a shorter RC decay time constant than the 
long cylindrical crystal of KDP.) 

The next two optical elements include a quarter wave plate and a horizontal polarizer. 
Except for its orientation, the polarizer is identical to that described above. The quarter wave 
plate is composed of quartz and has a diameter of 20 mm. As discussed in section 2, the quarter 
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wave plate changes the transmission profile T(V) from that given in (9) to that given in (10) 
when the optics are perfectly aligned and reflections are negligible. In effect, the quarter wave 
plate improves the electro-optic sensitivity T'(V) at a crystal voltage V = 0. This is important 
since, before a lightning event, the voltage on the crystal is very nearly zero (very small and 
random stray voltages undoubtably exist across the crystal at all times). A subsequent change in 
crystal voltage AV due to the lightning produces a change in optical transmission given by 
AT = T'(Q)AV for AV small. Hence, in order to get large (measurable) changes in transmission, 
it is desirable to have T'(0) large. As will be detailed in section 4, we have found that small 
misalignments in the optics produce sufficiently large values of T'( 0) when the quarter wave 
plate is not inserted. 

The final component shown in Figure 3a is the photodetector. It has an active detector area 
of 100 mm 2 and was connected directly to the input of the SR570 current preamplifier with a 
short (30.5 cm) 50 Q plenum coaxial cable. The amplifier (not shown) has a low end sensitivity 
of 1 volt output per 1 pA current input (written as 1 pA/volt). It also has valuable input offset 
adjust buttons that allowed us to accurately null the transmitted laser light (at V— 0) so we could 
increase amplifier gain without over ranging. We typically collected lightning data at either the 1 
//A/volt or 50 nA/volt gain setting, depending on storm range. The amplifier also offered 
considerable band pass filtering capabilities. Lightning data were collected using the 10 - 30 Hz 
band-pass filter, primarily to avoid dc drift and 60-cycle noise sources. Application of the band- 
pass filter contributed to shortening the overall, or “effective RC ” decay time constant of the 
sensor to a value of tens of milliseconds. The specific value of the effective RC decay depended 
on the particular amplifier gain setting that was chosen (see section 5 below for additional 
comments on sensor RC decay and specific examples). 
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Figure 3(b) shows the protective housing that the optical system is placed within. The 
housing is constructed of a 50.8 mm thick polypropylene material that affords very good thermal 
insulation from excessive heating due to typical summer conditions in Alabama. Ample room 
inside the box allowed for the insertion of cooling packs as an added measure to prevent excess 
heating to the laser and crystal; details of KDP refractive index changes with temperature are 
provided in [26]. Methods more sophisticated than applying cooling packs are the so-called 
“quenching conditions” discussed in [27] that remove unwanted birefringences in crystals 
affected by general changes in environmental conditions. Fortunately, under good thunderstorm 
conditions the weather was cloudy and not excessively hot, particularly during periods of 
precipitation. In these situations, we saw no discemable difference in sensor output with or 
without cooling packs. 

The polypropylene material had sufficient weight and shock absorption to minimize the 
effects of vibrations due to wind/precipitation/thunder that might disturb the optical alignment. 
Great care was spent assuring that the box was watertight; the box top lid overhangs the sides 
and a watertight O-ring lines the area where the lid rests. The outside of the box was coated with 
a durable water-repellent paint and the color was white to reflect sunlight (again preventing 
excessive heating inside the box). The inside of the box was coated with a flat black paint to 
minimize scattering of laser light that might be reflected from optical elements. 

On the outside top portion of the box lid is mounted a teflon rod that supports a circular 
aluminum antenna plate. A small bore hole was drilled through the teflon rod so that a teflon- 
coated wire could inconspicuously connect the antenna plate to the left ring electrode of the 
crystal. The other crystal ring electrode was grounded to the optical breadboard, as was the 
current preamplifier. The optical breadboard and aluminum ground plane (mounted on top of the 
box lid) were both connected to a standard 8 ft copper ground rod that was driven into the Earth 
a few inches from the housing. A small teflon “rain cap” was placed over the teflon support rod 
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before bolting the antenna plate down; this prevents streams of water from shorting the antenna 
to ground potential. 

In summary, the lightning-caused A E excites the sensor antenna and raises it to a transient 
voltage, AV, measured between the antenna plate (the left-side crystal electrode) and ground (the 
right-side crystal electrode); here we have neglected the small capacitance between the antenna 
plate and the left crystal electrode. The voltage AV across the crystal changes the crystal’s 
refractive index and polarization characteristics, and hence the transmission of the optical system 
changes by an amount AT. With laser power P () , the change in transmitted power is simply 
AP = ATP (I which is measured using a high gain photodetector/preamplifier system. Therefore, 
one formally has the chain AE -> AV -» AT -> AP so that values of AP can be directly related 
to the values of AE through calibration (see section 5 below). 

4. Alignment 

The expressions for the optical transmission, T, given in (9) and (10) assume that the optical 
elements are perfectly aligned. Since in this study the Pockels cell is used as a detector to sense 
the slightest antenna voltages, we have found that even small misalignments of the crystal and/or 
quarter wave are important effects. 

The theory of transmission through birefringent media is well developed and the 
transmission through a crystal at an oblique angle can be derived by exact or approximate means 
([19], [20], [28], [29]). When a collimated beam of light enters a birefringent medium at an 
oblique angle, the transmitted field consists of two beams that propagate in two distinct 
directions. The direction of a refracted beam depends on its angle of incidence and the refractive 
index, and the refractive index in turn depends on the angle of incidence. An exact solution to 
this problem is complicated, but a great simplification can be made by noting that there is very 
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little divergence in the transmitted beams, particularly for small angles of incidence. We will 
assume that the transmitted field consists essentially of a single beam whose direction is 
determined by the ordinary index of refraction. 

As defined in section 2, the crystal axes ( x , y, z) are parallel to the breadboard axes (X, Y, Z) 
when the optics are perfectly aligned. Suppose that the crystal is misaligned by an angle —0 

about an axis h = -sin^X + cos^Y, where <f> is measured from the positive A'-axis. Inside the 
crystal, both refracted rays are assumed to make an angle 0with the crystal’s extraordinary axis, 
where sin© = n 0 sin (9. A formula for retardation may be obtained by examining the index 
ellipsoid. 

The equation for the index ellipsoid and expressions for related variables can be 
summarized as 

-2 c A \ 

n„ B 0 

r r -N-r = l, r = xx + >y + zz, N= S n~ 2 0 , S — e z r 63 . (11) 

V o 0 n' 1 ) 

Ideally, light should pass along z; however, in the misaligned case the optical path is 
along sin#(cos^x + sin^y) + cos#z. Such a direction may be obtained by rotating z about h 
through an angle 0. This may be done with a rotation matrix whose components are obtained 
from the dyadic 

M = [I - (1 - cos 6>)(I - hh) + sin 6>h x i], (12) 

where f is the identity dyadic. In a rotated system, the equation of the index ellipsoid and 
related variables are given by 

r' 7 -N'-r' = 1, M r' = r, N' = M T -N-M, (13) 

and the optical path becomes z r . 
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Polarizations and indices of refraction may be obtained from the components of N 2 given 


by: 


N 


2 — 



0 ] . 

y) ’ 


a = ( A 2 + A 2 )h„“ 2 + A 2 n/ 2 + 2 A X A 2 £ 

/3=A 2 (A l + A 4 )n 0 ~ 2 +A 3 A 5 n e ~ 2 + (A ] A 4 +A 2 )s 
r=(Aj + A] )n 2 + A]n 2 + 2A 2 A 4 £ ; 

A, = cos#cos 2 ^ + sin 2 ^, A 2 = (cos#-l)cos^sin^, A 3 = sin#cos^, 
4, =cos#sin 2 ^ + cos 2 ^, ,4 5 = sin#sin^. 


(14) 


where the intersection to the plane normal to the beam direction and the index ellipsoid is the 


ellipse r 2 7 -N 2 -r 2 ' =1 and r 2 =x'x' +y'y ' . The refractive indices and allowed polarizations 
may be obtained from the eigenvalues and eigenvectors of this matrix. Performing these 
calculations yields the two eigenvalues 

X ± = y {2 n~ 2 + sin 2 OTt - sin 2 0 sin 2<f>S 

± [sin 4 Q7t +2(l + cos 2 #)sin 2 #sin2 <j>1t S 

. (15) 

+ (sin 4 #sin 2 2<p + 4cos 2 0)S 2 J j 

««„‘ 2 + 0 2 {^t - sin2^5)/2 ± [# 4 % 2 /4 + # 2 sin2 </>7t £+£ 2 f 
for small 6, where % = n~ 2 - n ~ 2 . The retardation is T = 2nd sec# A/ X, where 

A = n x -n f = l/VT - \/4k ~ T j4Vd 2 sin2^E + E 2 . (16) 

Eigenvectors are given by 

v'_ = cos^x' -sin#y', v' + = sin<Jx' + cos^y', (17) 


where tan25 =2fi/(y -a) . The + and - signs indicate the fast and slow axes, respectively. 
Consistent with the approximations used, the field exiting the second polarizer has the form 
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( 18 ) 


F 0 [(V + • X)(v + • Y)cos<y r + (v. • x)(v_ • y)cos(6^ - T)]x, 

V_ = cos£X-sinc>Y, V + = sin£X + cos£Y. 

Obtaining the transmission as for (9), we have 

r = sin 2 2<Wr/2, (19) 

where sin 20 =2/?/ ^/(a-/) 2 + 4/? 2 . If a quarter wave plate is used, similar reasoning leads to 
the formula 

T = \[\ + sin20 sinT], (20) 

The choice of sign for the square root corresponds to the quarter wave plate installed in such a 
way as to increase transmission for increasing voltage (e z being taken to be positive in the 
positive direction of the Z-axis). 

In order to provide for a most general misalignment of the crystal, we allow for a rotation of 
the crystal about Z by an angle r. Results may be obtained by rotating the orientations of the 
polarizers and quarter wave plate by — t. Transmissions in this case are given by (19) and (20) 
with 8 replaced with 8 + t. 

Using the expressions in (19) and (20), we have computed the theoretical transmission of an 
aligned and misaligned optical system with or without a quarter wave plate. The results are 
provided in Figure 4. We have estimated the small misalignments expected under normal 

laboratory conditions, i.e., 0 = 1.4° so that 0 = sin -1 (1.5074 sin 1.4) = 2.1 1° about an axis h 

oriented with angle (p = 89.75° . 

Note that when there is no quarter wave plate and the optics are aligned (19) reduces to (9) 
and a simple sine squared dependence results (solid line. Figure 4(a)). When the optics are 
misaligned, the net transmission increases slightly, but more importantly, the sensitivity T'(V) 
increases substantially at and near V = 0 (dashed line, Figure 4(a)). When the quarter wave is 
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inserted, the transmission increases for both aligned (dotted line, Figure 4(b)) and misaligned 
(dot-dashed line, Figure 4(b)), but the misalignment no longer improves sensor sensitivity at and 
near V= 0, or at any larger values of V. 

To further clarify the variation of sensor sensitivity for small changes in crystal voltage, we 
have numerically computed the change in transmission AT = T(V)-T(0 ) for a small crystal 
voltage change AV = V-0= V all as a function of V. The functional dependence AT = f(V) 
is provided in Figure 5 for aligned and misaligned optics with or without a quarter wave plate. 
We see that the addition of the quarter wave plate necessarily improves sensor sensitivity no 
matter if the system is aligned or misaligned. However, because the quarter wave plate produces 
additional reflective losses to the optical system, because its inclusion produces additional noise 
sources during high-gain operation of the lightning detector, and because the slope of the lines in 
Figure 5(b) have magnitudes that are not vastly different, we did not include the use of the 
quarter wave plate during normal calibration and operation of the lightning sensor. 

5. Calibration 

The calibration of the lightning detector was straightforward. On the top lid of the optical 
housing shown in Figure 3a, four threaded holes were drilled (one at each corner of the square 
lid). A nylon rod was screwed into each hole in order to support a large circular aluminum 
calibration plate. An Avtech high voltage wave function generator (HVWFG) was connected 
between the calibration plate and the lid-mounted ground plane. This provided a uniform vertical 
electric field change for the lightning detector antenna plate located between the calibration plate 
and ground plane. We also used a standard (i.e., lower voltage) BK-Precision WFG to check 
low-end sensor detectability. 
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The risetime of a lightning return stroke field change is typically 60 //s. The Avtech 
HVWFG could supply square waves having amplitudes of ±500 volts with frequencies ranging 
from near dc to 3 kHz and with a peak-to-peak risetime < 35 /is. The BK Precision WFG could 
supply square waves having amplitudes of ±12 volts with frequencies ranging from near dc to 
several megahertz, and a risetime of about 2 ms when operated at 0.5 Hz. In order to get an 
accurate calibration, it is important to have waveforms that rise to peak in milliseconds or less 
because of the characteristic RC decay of the crystal. In effect, the vertical displacement current 
density e o dEjdt produced by the calibration apparatus or a lightning event must be large 
enough to place a significant amount of charge onto the left crystal electrode before leakage 
currents in the crystal drain the crystal voltage to ground. The electric field line distribution 
inside an electro-optic parallelepiped crystal is provided in [30] and evidence for field-driven 
leakage currents in electro-optic materials is given in [3] and [24]. 

Because we have picked a crystal that has an RC decay of ~1 s, the leakage current to 
ground is slow. This allows the crystal to acquire a larger voltage, and hold a voltage for a longer 
period of time. Furthermore, since electro-optic crystals can be modulated into the gigahertz 
region, they respond quickly to changes in voltage. Hence, the crystal responds quickly to 
lightning transients, acquires a sufficient voltage, and holds the voltage for a time that is long in 
comparison to the (tens of milliseconds) decay time of the photodetector/amplifier system. The 
amplifier effectively samples the crystal voltage substantially faster than the crystal voltage has a 
chance to decay, thereby optimizing overall sensor sensitivity. 

If it were not for our desire to remove 60-cycle noise sources, we could replace the applied 
band-pass filter with a high-pass filter (> 10 Hz) to improve overall sensor bandwidth. A narrow 
filter for removing 60-cycle noise in conjunction with the high-pass filter would clearly be 
optimum, but is not possible with the SR570 electronics alone. 
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Figure 6 shows some specific examples of sensor output in response to both high and low 
voltage WFG-simulated lightning electric field changes; the amplifier gain was set to 1 //A/volt 
and the 10 - 30 Hz band-pass filter was used. The train of sensor output pulses provided in 
Figure 6(a) are due to a ±2083 V/m square wave calibration input; this field change magnitude is 
typical of lightning that is a few kilometers away from the observer. Note that the output 
waveforms decay back to zero with an e-folding time (effective RC decay) of ~14 ms. The 
amplifier then “overshoots” until it recedes back completely to zero about 1 54 ms later. At a gain 
of 50 nA/volt, there was no overshooting and the e-folding time was ~27 ms. 

Figure 6(b) shows the output for a square wave input of only ±30.6 V/m. The low-end 
sensitivity and signal-to-noise shown in Figure 6(b) can be further improved by further 
optimizing the misalignment of the optical system. Additional details related to the general 
accuracy and limits of precision of electro-optic voltage sensors is provided in the U.S. National 
Bureau of Standards (NBS) report given in [31]. 

We performed these type of calibration tests for a variety of simulated lightning A E values 
and amplifier gains when the sensor was deployed outside of building 4614 at the NASA MSFC. 
Figure 7 shows a sample linear response calibration plot acquired over a large range of A E values 
when the gain was set to 1 //A/volt. Because the KDP crystal half wave voltage V n is about 9000 
volts, the lightning detector’s full dynamic range well exceeds that shown in Figure 7. In 
addition, because of the misalignment of optical elements, the value of S in (19) differs from 
tz/4 so that we are able to discriminate between positive and negative lightning A£s, i.e., the 
linear calibration curve extends to negative values of A E. 

6. Storm Data 


20 



The summer of 1998 in northern Alabama was marked by substantially dry periods that led 
to eventual fire outbreaks in September. Nonetheless, we were able to deploy the lightning 
detector on five separate occasions for storms that occurred on June 30, and July 3, 20, 21, and 
24. The storm occurring on July 20 resulted in a lightning death to a citizen of northern 
Huntsville, Alabama that occured about 45 minutes prior to the start of our data collection 
period. 

The National Lightning Detection Network (NLDN) operated by Global Atmospherics Inc. 
(GAI) located in Tucson, Arizona detected widespread cloud-to-ground activity on June 30 over 
much of northern Mississippi, Alabama, and Georgia, and over most of Tennessee. Several 
cloud-to-ground lightning events that occurred within a few kilometers of our deployment site 
were identified by eye and were time-correlated with transient sensor outputs that we displayed 
in real time on the Pentium computer. We used Microcal Origin 4.0 data acquisition software in 
“scope” mode to look at the digitized sensor output. We also viewed the direct analog output of 
the sensor on a conventional 40 MHz oscilliscope. Our sensor detects both ground and cloud 
discharges, but for these daytime thunderstorms it was more practical to time-correlate events 
with ground discharges, i.e., it was easier to see the well-defined channel of a ground stroke 
rather than the faint diffuse emission from cloud-obscured discharges. 

Figure 8(a) shows the output from the electro-optic lightning detector for an active 2-min 
period during the June 30 storm (the plot abscissa is in seconds after 4:48:30.5 PM CST, i.e., the 
time the sensor was turned on). These data were collected using a gain of 50 nA/volt, or a factor 
of 20 higher than that shown in the calibration curve in Figure 7. By looking southward from our 
deployment site during this storm, we were able to visually confirmed over 10 lightning events 
that gave coincident field excursions in our recorded sensor output. Two small amplitude 
lightning flashes are indicated with pointers in Figure 8(a), and a more time-resolved plot of each 
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flash is provided in Figure 8(b) and 8(c). These particular flashes were confirmed by eye as 
cloud-to-ground events. Several rapid variations are evident in each flash that we believe are 
individual strokes. The duration between the peaks of these excursions is less than -100 ms 
which is characteristic of the interstroke time interval shown in Figure 3 of [2], There may even 
be evidence of a small amplitude “K” change just after 3160.4 s in Figure 8(c). There are some 
positive excursions starting at -3122.7 s in Figure 8(b) that might be associated with stepped 
leader activity or other local atmospheric electrical effects. We believe that the negative 
excursion before the series of positive strokes is associated with an upward in-cloud moment 
change; a similar negative excursion also occurs in the flash shown in Figure 8(c). 

As we have stated, vibrations due to wind, precipitation, and thunder can produce noise in 
the sensor output and these affects are evident for small amplitude discharges. Figure 9 shows 
some time-resolved, single pulse, large amplitude events that occurred earlier during the June 30 
storm. There are several low amplitude reverse-polarity field changes associated with distant 
lightning events, and two larger amplitude events that are probably much closer to the sensor. It 
is possible that the reverse-polarity events are associated with cloud discharges that occur 
beyond the “reversal distance”, i.e., the surface field verses range from a tripolar thundercloud 
charge distribution reverses sign at the “reversal distance”, or about 10 km [1]. 

The saturation output voltage of the SR570 amplifier is ±5 volts which was occasionally 
exceeded by some lightning events. In some cases, we were forced to decrease sensor gain to 
avoid frequent pulse saturation in the data stream. An RS232 cable connection between the 
amplifier and our indoor computer allowed us to remotely adjust sensor gain as the storm cells 
moved or developed relative to the deployment site. 
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7. Summary 


In the future, we would like to further experiment with and optimize the electro-optic 
lightning sensor and deploy it again to collect more storm data. As we have stated, the 
dimensions of our sensor could be scaled down significantly. The laser and all optical elements 
could be precision-mounted within a single sealed tube only a few inches long and wide. 
Alternately, a small crystal attached to polarization maintaining (PM) fibers and an antenna 
could be the sole components deployed into the thunderstorm environment; the photodetector, 
current amplifier, laser, and remaining optics could all reside inside a building along side the 
data acquisition computer. Or, crossed polarizers could be mounted to the crystal and regular 
fiber optic cables used. In any case, a single unit reduced size system would be less susceptible 
to vibration-induced noise sources. 

Size reduction, as well as improvements in low end sensitivity and signal-to-noise are 
fundamental to improving the attractiveness of electro-optic applications in the discipline of 
lightning research. It is possible to arrange optics in such a way as to pass the laser light multiple 
times through the crystal to increase the total birefringence and thereby substantially improve 
sensor gain. We are also interested in inspecting other crystal types that offer certain advantages 
such as bismuth germanate (Bi4Ge30]2) described in [32] and Hexamine ((CH 2 ) 6 N 4 ) described 
in [33], 

Ultimately, the electro-optic technology needs to be placed beside a conventional electric 
field mill such as the type operated by the NASA Kennedy Space Center (KSC) and the U.S. Air 
Force Eastern Range (ER) located in east central Florida. This would allow one to further inspect 
the sensitivity, accuracy, signal, and noise characteristics of the electro-optic sensor under a 
variety of environmental conditions and lightning source stimuli. 
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In the future, lightning field changes derived from a ground-based network of electro-optic 
sensors that cover a region of a few hundred square kilometers could be used to retrieve the 
charges deposited by a lightning. This has already been achieved using the KSC-ER field mill 
network, but with the added difficulty and errors associated with detecting and computing values 
of A E from the continuous £-data ([33] - [37]). 
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Figure Captions 


Figure 1. Typical electro-optic Pockel cell light amplitude modulator [adapted from Yariv and 
Yeh, 1984], 

Figure 2. Overview of the electro-optic lightning detector system; A E denotes a thunderstorm 
electric field change produced by a cloud-to-ground lightning event shown to the left. 

Figure 3. The basic components of the electro-optic lighting detector: (a) optics, and (b) the 
optical housing and antenna. 

Figure 4. Comparison of the transmission T for aligned and misaligned optics (a) without a 
quarter wave plate and (b) with a quarter wave plate. 

Figure 5. The change in transmission AT for a voltage change AV = V - 0 = V across the crystal 
all as a function of V, with and without a quarter wave plate for (a) aligned optics and (b) 
misaligned optics. 

Figure 6. Some calibration results for (a) a typical lightning field change of 2083 V/m and (b) a 
very small 30.6 V/m field change that would be produced by a distant lightning source. The 
band-pass filter of the amplifier is 10 - 30 Hz and the amplifier gain is 1 //A/volt . 
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Figure 7. Sample calibration curve across a wide dynamic range of A£ values. A linear 
regression fit to the data points shown is A E = -21.130 + 371.495AF; the linear regression 
constant is r = 0.99902. The band-pass filter of the amplifier is 10 - 30 Hz and the amplifier gain 
is 1 //A/volt. 

Figure 8. Example of lightning data results for a storm that occurred on June 30, 1998: (a) for a 
2-min period, (b) for a cloud-to-ground lightning flash within the 2-min period, and (c) another 
cloud-to-ground flash in the 2 minute period. The band-pass filter of the amplifier is 10 - 30 Hz 
and the amplifier gain is 50 nA/volt, or a factor of 20 higher than that used in the calibration plot 
of Figure 7. 

Figure 9. Additional examples of larger amplitude time-resolved lightning field changes 
collected from the June 30, 1998 storm. Several small field excursions exist due to distant 
lightning. The band-pass filter of the amplifier is 10 - 30 Hz and the amplifier gain is 50 nA/volt, 
or a factor of 20 higher than that used in the calibration plot of Figure 7. 
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